Moving Vortex Phases, Dynamical Symmetry Breaking, and Jamming for Vortices in 

Honeycomb Pinning Arrays 
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We show using numerical simulations that vortices in honeycomb pinning arrays can exhibit a 
remarkable variety of dynamical phases that are distinct from those found for triangular and square 
pinning arrays. In the honeycomb arrays, it is possible for the interstitial vortices to form dimer 
or higher n-mer states which have an additional orientational degree of freedom that can lead to 
the formation of vortex molecular crystals. For filling fractions where dimer states appear, a novel 
dynamical symmetry breaking can occur when the dimers flow in one of two possible alignment 
directions. This leads to transport in the direction transverse to the applied drive. We show that 
dimerization produces distinct types of moving phases which depend on the direction of the driving 
force with respect to the pinning lattice symmetry. When the dimers are driven along certain 
directions, a reorientation of the dimers can produce a jamming phenomenon which results in a 
strong enhancement in the critical depinning force. The jamming can also cause unusual effects 
such as an increase in the critical depinning force when the size of the pinning sites is reduced. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 



Vortex matter in type-II superconductors has been ex- 
tensively studied as a unique system of many interacting 
particles in which noncquilibrium phase transitions can 
be accessed readily [ll, 1^ 111 S IE 0, i- In the ab- 
sence of driving or quenched disorder, the vortex-vortex 
interactions favor a triangular crystalline ordering. If the 
sample contains sufficiently strong quenched disorder in 
the form of randomly placed pinning sites, the vortex lat- 
tice ordering can be lost as the vortices adjust their posi- 
tions to accommodate to the pinning landscape [U, [3) 9 ■ 
Under an applied drive such as the Lorentz force from 
a current, the vortices remain immobile or pinned for 
low driving forces; however, there is a threshold applied 
force above which the vortices begin to move over the 
quenched disorder. For strong disorder, the initial mov- 
ing state is highly inhomogeneous with the vortices flow- 
ing in meandering and fluctuating channels, and there is 
a coexistence between pinned vortices and flowing vor- 
tices [J, 0|- At higher drives the vortices move more 
rapidly, the effectiveness of the quenched disorder is re- 
duced, and the fluctuations experienced by the vortices 
become anisotropic due to the directionality of the ex- 
ternal drive The vortex- vortex interactions become 
more important at the higher drives when the quenched 
disorder becomes ineffective, and a dynamical transition 
can occur into a moving smectic state where the vortices 
regain partial order in one direction 0, @, 0] • Here, the 
system has crystalline order in the direction transverse to 
the vortex motion and liquid-like order in the direction of 
vortex motion. Depending on the dimensionality and the 
strength of the pinning, it is also possible for the vortices 
to reorganize in both directions at high drives to form 
a moving anisotropic crystal 0, H, 0, H, Q . The ex- 
istence of these different phases and transitions between 
the phases can be inferred from signatures in transport 



[1] and noise fluctuations [13, [Tlj , and the moving phases 
have also been imaged directly using various techniques 

In addition to the naturally occurring randomly placed 
pinning sites, it is also possible for artificial pinning sites 
to be created in a periodic structure [l3|. Recent ad- 
vances in nanostructuring permit the creation of a wide 
variety of periodic pinning landscapes where the period- 
icity, shape, size and density of the pinning sites can be 
well controlled. Distinct types of pinning arrays such as 
square [H [lljlljll, [13, 111 [H, [M triangular [2T, ^ , 
rectangular [23l |24| . honeycomb [25l [26j . kagome {2^ , 
quasicrystalline [IJ , and partially ordered structures 
have been created. In these arrays the type of vortex 
structure that forms is determined by whether the vor- 
tex lattice is commensurate with the underlying pinning 
array. Commensurate arrangements appear at integer 
multiples of the matching field B^, which is the mag- 
netic field at which the vortex density matches the pin- 
ning density, and in general, ordered vortex states oc- 
cur at matching or rational fractional values of B/B^ 
[H, [11, [H, [H, [33, [U, HI. In samples where only one 
vortex can be captured by each pinning site, the vor- 
tices that appear above the first matching field sit in the 
interstitial regions between the pinning sites, and these 
interstitial vortices can adopt a variet y of crystalline con- 
figurations [11, m, [H, [13, IS [M mM m ■ 

Since a number of distinct ordered and partially or- 
dered vortex states can be created in periodic pinning 
arrays, a much richer variety of dynamical vortex behav- 
iors occur f or p eriodic pinning than for random pinning 

II' III ssMJii, m, [13, m, m, no, m 

44 I45I l46l I47I l48l l49j |. Several of the dynami- 
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cal phases occur due to the existence of highly mobile 
interstitial vortices which channel between the pinned 

vortices [H, [H IH, [H, HI M M M \M, M 

As a function of applied drive, various types of mov- 



2 



ooooooooooo 
oooooooooooo 

ooooooooooo 
oooooooooooo 

ooooooooooo 
y oooooooooooo 

ooooooooooo 
oooooooooooo 

ooooooooooo 
oooooooooooo 

ooooooooooo 
oooooooooooo 

ooooooooooo 



(a) 



o 


o 


o 


o 


o 


o 


o 







o 


o 


o 


o 





o 








o 








o 





o 







o 


o 


o 


o 





o 





o 


o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 







o 





o 


o 





o 








o 








o 





o 







o 


o 


o 


o 





o 





o 


o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 





(b) 



(c) 



(d) 



FIG. 1: (a) Pinning site locations (open circles) for a triangu- 
lar pinning array, (b) Pinning site locations for a honeycomb 
pinning array constructed from the triangular array in (a) by 
removing 1/3 of the pinning sites, (c) The pinning site lo- 
cations and vortex positions (dots) for a honeycomb pinning 
array at B/B^ = 1.5. The overall vortex lattice order is tri- 
angular, (d) The pinning site locations and vortex positions 
for a honeycomb pinning array at B/B^ — 2.0, where two 
vortices are captured at the large interstitial sites and the re- 
sulting dimers all have the same orientation. Here Fp — 0.85, 
Rp = 0.35A, and for the honeycomb array Up = 0.3125/A^. 



ing phases occur, including interstitial vortices moving 
coherently between the pin ning sites in one-dimensional 
paths m, H M,^,^,^ or periodically mod- 

ulated winding paths [33, 35, 44, 45], disordered reg imes 
where the vortex motion is liquidlike [H, |3^ l39l . |4^ , 
and reg imes where vortices flow along the pinning rows 
[H, [li, m, 113, nil . Other dynamical effects, such as 
rectification of mixtures of pinned and interstitial vor- 
tices, can be realized when the periodic pinning arrays 
are asymmetric (49j . 

Most of the studies of vortex ordering and dynamics in 
periodic pinning arrays have been performed for square 
and triangular arrays. Experiments with honeycomb and 
kagome pinning arrays revealed interesting anomalies in 
the critical current at nonmatching fields which are as 
pronounced as the anomalies observed at matching fields 
in triangular pinning arrays [25l . [26j . A honeycomb pin- 
ning array is constructed by removing every third pin- 
ning site from a triangular pinning array, producing a 
periodic arrangement of triangular interstitial sites. In 
Figs. l(a,b) we illustrate a triangular pinning array along 
with the honeycomb pinning array that results after the 
removal of one third of the pinning sites. The match- 



ing anomalies in the experiments coincide with fields 
B/Bfp = m/2, with m an integer. At these fields, the 
vortex density would match with the regular triangular 
pinning array. At the matching anomalies for m > 2, a 
portion of the vortices are located in the large intersti- 
tial regions of the honeycomb lattice [l^, as illustrated 
in Fig. 1(c) for B/B^ = 1.5. The overall vortex lattice 
structure is triangular and a strong peak in the depinning 
force occurs at this field (soj . 

Recently, we used numerical simulations to demon- 
strate that vortices in honeycomb pinning arrays have 
a rich equilibrium phase diagram as a function of vor- 
tex density ^50] , with matching anomalies at integer and 
half integer matching fields that are in agreement with 
experiments. The large interstitial sites created by the 
missing pinning sites can capture multiple interstitial 
vortices which form cluster states of n vortices. For 
1.5 < B/Bfj, < 2.5, dimer states with n = 2 form, while 
for higher fields trimer and higher order n-mer states 
form. At the integer and half-integer matching fields, 
the n-mer states can assume a global orientational or- 
dering which may be of ferromagnetic or antiferromag- 
netic type; herringbone structures can also form, similar 
to those observed for colloidal particles on periodic sub- 
strates |5L 52, 53, 54, 55] and molecules on atomic sub- 
strates |56| . These orientationally ordered states have 
been termed vortex molecular crystals. Certain vortex 
molecular crystals have ground states that are doubly or 
triply degenerate, such as the dimer state illustrated in 
Fig. 1(d) at B/B^ = 2.0 where the dimers align in one 
of three equivalent directions [13] ■ As the temperature 
is increased, the n-mers undergo a transition from an or- 
dered state to an orientationally disordered state in which 
the n-mers rotate randomly but remain confined to the 
interstitial pinning sites. The rotating states have been 
termed vortex plastic crystals. At matching fields where 
vortex plastic crystals form, the anomalies in the critical 
current disappear [50| . The predictions from the simu- 
lations are in general agreement with the experimental 
observation of the loss of certain higher order matching 
anomalies at higher temperatures [26|. The formation of 
n-mers that can be aligned along degenerate symmetry 
directions has also been predicted for kagome pinning ar- 
rays where every other pinning site is removed from every 
other row of a triangular lattice (sol . [STj . 

The formation of dimer states in the honeycomb pin- 
ning array leads to a variety of novel dynamical phases, 
including a spontaneous dynamical symmetry breaking 
effect in which the moving vortices organize into one of 
two equivalent states which have a component of trans- 
lation perpendicular to the applied drive in either the 
positive or negative direction [58|. The transverse re- 
sponse appears when the external driving force is applied 
halfway between the two symmetric directions of aligned 
dimer motion. The dynamical symmetry breaking oc- 
curs when the equilibrium ground states have no global 
symmetry breaking. At B /B^ — 2.0, the ground state is 
symmetry broken and the dynamical moving state has 
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the same broken symmetry as the ground state. For 
incommensurate fiUings, when the dimer ahgnment is 
disrupted, there is no global symmetry breaking in the 
ground state, and instead a dynamical symmetry break- 
ing occurs due to the applied drive. 

In this work we map the dynamical phase diagram for 
vortices in honeycomb arrays. We focus on the states 
1.5 < B/B^ < 2.5 to understand where dynamical sym- 
metry breaking occurs and to examine what other types 
of moving phases are possible. We study how the dy- 
namical phases change for driving along different axes 
of the pinning lattice. We find that very different kinds 
of dynamics occur when the driving direction is varied, 
and that the value of the depinning threshold is strongly 
directionally dependent. We also find that a novel jam- 
ming phenomenon can occur due to the formation of the 
dimer states. For certain directions of drive, the dimers 
are anti-aligned with the drive, causing the dimers to be- 
come blocked in the interstitial regions. 

Although our results are specifically for vortices in 
type-II superconductors, the general features of this work 
should also be relevant for other interacting particle sys- 
tems where a periodic substrate is present. Examples of 
such systems include vortices on periodic substrates in 
Bose-Einstein condensates (BEC), where different kinds 
of crystalline phases can occur which depend on the 
strength of the substrate [H, It should be pos- 

sible to observe different types of vortex flow states in 
BEC systems [6l!|. Our results are also relevant for col- 
loids on periodic substrates, where an orientational or- 
dering of colloidal molecular crystals occurs which is very 
similar to that of the vortex molecular crystal states 
[m m m HI HE [11. other related systems include 
charged balls on periodic substrates [63*1 and models of 
sliding friction 64\. 

II. SIMULATION 

We use the same simulation employed in the previous 
study of vortex equilibrium states in honeycomb pinning 
arrays [131 . We consider a 2D system of dimensions = 
L and Ly = L with periodic boundary conditions in the x 
and y directions. The sample contains Ny vortices, giving 
a vortex density of = N^/L"^ which is proportional to 
the external magnetic field. In addition, there are Np 
pinning sites placed in a honeycomb arrangement with a 
pinning density of Up = Np/L"^ . The field at which the 
number of vortices equals the number of pinning sites is 
defined to be the matching field B^. 

The dynamics of vortex i located at position is gov- 
erned by the following overdamped equation of motion: 

^^ = Fr+Fr+F,,+Ff. (1) 

Here the damping constant is 77 = 4%d/2'KS^^ p[^, where 
d is the thickness of the superconducting sample, ^ is 
the superconducting coherence length, pat is the normal 



state resistivity of the material, and 0o = h/2e is the 
elementary flux quantum. The vortex-vortex interaction 
force is 

Fr = E/oi^i(^)R., (2) 

where Ki is the modified Bessel function, A is the London 
penetration depth, /g — (I)q/{2tt^oX^), Rij = |R.i — Rj| is 
the distance between vortex i and vortex j, and the unit 
vector Rij = (R— Rj)/i?y . In this work all length scales 
are measured in units of A and forces in units of fo- The 
vortex vortex interaction decreases sufficiently rapidly at 
large distances that a long range cutoff is placed on the 
interaction force at Rij — 6A to permit more efficient 
computation times. We have found that the cutoff does 
not affect the results for the fields and forces we consider 
here. 

The pinning force F^^ originates from individual 
nonoverlapping attractive parabolic traps of radius Rp 
which have a maximum strength of Fp. In this work we 
consider the limit where only one vortex can be captured 
per pinning site, with the majority of the results obtained 
for Rp — 0.35A. The exact form of the pinning force is: 

Fr = - E fo ) R^.'^ (^^^Sr^) 

Here, i?-^"* — |Ri — R^^''|, R''^ ^ is the location of pinning 

site k, the unit vector R^^^^ = (R^ - R[^V^lfc\ and 9 is 
the Heaviside step function. 

The external drive = Fd/qFd represents the 
Lorentz force from an applied current J x B which is 
perpendicular to the driving force and is applied uni- 
formly to all the vortices. We apply the drive at various 
angles to the symmetry axes of the honeycomb pinning 
array. The thermal force Ff originates from Langevin 
kicks with the properties (Ff ) = and (Ff (t)Fj(i')) = 
2rjkBT5ij5{t — t'). Unless otherwise noted, the ther- 
mal force is set to zero. The initial vortex configura- 
tions are obtained by simulated annealing, and the ex- 
ternal force is then applied gradually in increments of 
AFd = 0.0002 every 1000 simulation time steps. For 
the range of pinning forces used in this work, we find 
that this force ramp rate is sufficiently slow that tran- 
sients in the vortex dynamics do not affect the overall 
velocity-force curves. We obtain the velocity- force curves 
by summing the velocities in the x (longitudinal) direc- 
tion, (V^) = X]i=i ' and the y (transverse) di- 
rection, {Vy) = N-'^Y.t • y, where = dKi/dt. 
In Fig. l(c,d) we illustrate the pinning sites and vortex 
configurations after simulated annealing for B/B^ = 1.5 
[Fig. 1(c)] and 2.0 [Fig. 1(d)]. Here = Ly = 24A and 
Up = 0.3125/A-^. In our previous work, Ref. [Ill, the 
drive was applied along the x-direction for the geometry 
in Fig. 1. 
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FIG. 2: (a) The average velocity in the a;-direction (Vx) vs 
external driving force Fd for the honeycomb pinning array 
from Fig. 1(d) at B/B4, = 2.0 with Fd = Fdx. (b) The cor- 
responding average velocity in the {/-direction {Vy} vs Fd- We 
observe four phases: the initial pinned phase (P), a symmetry 
broken phase (SB), a random phase (R), and a moving locked 
phase (ML). 



III. DYNAMICS AND TRANSVERSE 
RESPONSE FOR DRIVING IN THE 
LONGITUDINAL DIRECTION 

We first consider the case for driving in the x or lon- 
gitudinal direction, Fd = F^x., for the system shown in 
Fig. 1(d) with B/B^ = 2.0, Rp = 0.35A, and Fp = 0.85. 
In Figs. 2(a,b) we plot (I4) and (Vy) versus Fjj. At this 
filling there are four distinct dynamical phases, with the 
pinned (P) phase occurring at low Fjj. The depinning 
threshold Fc occurs near Fd = 0.14 when the interstitial 
vortices become depinned. For a system with random 
pinning and Fd = -Fdx, there would be no transverse 
velocity response; the system would have {Vy) = and 
only (Vx) would be finite. In contrast, for the honeycomb 
pinning array there is a finite velocity both in the posi- 
tive X direction and in either the +y or —y direction. In 
Fig. 2(b) the transverse response {Vy) is negative, indi- 
cating that the vortices are moving at a negative angle to 
the X axis for 0.14 < Fd < 0.37. Figure 3(a) illustrates 
the vortex motion at Fd — 0.25, where the vortices flow 
in one-dimensional paths oriented at —30° to the x axis. 
In Fig. 3(b) a snapshot of the vortex positions shows that 
the vortex lattice remains ordered in the moving phase, 
indicating that the vortices are flowing in a coherent man- 
ner. We term the phase shown in Fig. 3(a) the symmetry 
broken (SB) phase, since the flow can be tilted in either 
the positive or negative y-direction. 

At B/Bfj, = 2.0 and Fd = 0, the interstitial vortices 




W X (c) X (e) 



(H) 5 W) 5 i 

FIG. 3: The dynamics of the three moving phases from 
Fig. 2 for the honeycomb pinning array at B/B^ = 2.0 with 
F_D ~ Fd^- The vortex positions (filled circles), pinning site 
locations (open circles), and vortex trajectories (black lines) 
are shown in an 18A x 18A portion of the sample, (a) In the 
symmetry broken (SB) phase at Fd ~ 0.25, the interstitial 
vortices move along a —30° angle to the a;-axis while the vor- 
tices at the pinning sites remain immobile, (b) Vortex posi- 
tions only in the SB phase at Fd = 0.25, showing the ordering 
present in the vortex lattice structure, (c) In the random (R) 
phase at Fd = 0.42, the vortex motion is highly disordered 
with vortices pinning and repinning at random, (d) Vortex 
positions only in the R phase at Fd ~ 0.42 indicate that the 
vortex lattice is disordered, (e) In the moving locked (ML) 
phase at Fd ~ 0.65, all the vortices channel along the pinning 
sites, (f) Vortex positions only in the ML phase at Fd ~ 0.65 
reveal an anisotropic vortex lattice structure with different 
numbers of vortices in each row. 



form an aligned dimer configuration with a three-fold de- 
generate ground state in which the dimers can be oriented 
along the y-direction, as in Fig. 1(d), or along -1-30° or 
—30° to the a;-direction, as shown in previous work [ssj . 
When a driving force is apphed to the +30° or —30° 
ground states, the vortices depin and flow along -1-30° or 
—30°, respectively. In these cases, the symmetry break- 
ing in the moving state is not dynamical in nature but 
reflects the symmetry breaking within the ground state. 
If the dimers are initially aligned along the y-direction 
in the ground state, an applied drive induces an instabil- 
ity in the pinned phase and causes the dimers to rotate 
into the -1-30° or —30° directions, as we discuss in fur- 
ther detail below. In this case the symmetry breaking is 
dynamical in origin. 

In Fig. 2(a,b) we find pronounced oscillations in both 
(Vx) and {Vy) just above the depinning threshold Fc — 
0.14. These oscillations arc not intrinsic features but are 
due to the fact that at B/B^ = 2.0 the interstitial vor- 
tex lattice is perfectly ordered, so the interstitial vor- 
tices move in a coherent fashion as shown in Fig. 3(a). 
At depinning, the interstitial vortices are slowly mov- 
ing through a periodic potential created by vortices that 
remain trapped at the pinning sites. This periodic po- 
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FIG. 4: Time traces of vortex velocity at fixed Fd- Upper 
curves: Vx{t)\ lower curves: Vy{i). (a) The symmetry broken 
(SB) phase at Fd = 0.25 from Fig. 3(a,b). Here pronounced 
oscillations occur in both 14 and Vy as the vortices move in 
a coherent fashion, (b) The random (R) phase at Fd ~ 0.42 
from Fig. 3(c,d). In this case the transverse motion is lost 
and (Vy) = 0. Additionally, there are no correlated oscilla- 
tions, (c) The moving locked (ML) phase at Fd = 0.65 from 
Fig. 3(e,f). Vx has been shifted down for clarity. There is a 
weak oscillation in Vx due to the periodic substrate. Since the 
flow is strictly one-dimensional, as shown in Fig. 3(e), there 
are no fluctuations in Vy. 



tential causes the moving interstitial vortices to develop 
an oscillating velocity. In Fig. 4(a), the instantaneous 
time traces of the vortex velocity and Vy at constant 
= 0.25 show strong velocity oscillations. The os- 
cillations are also visible in Fig. 2 at low drives due to 
our choice of averaging time spent at each value of the 
driving current. If the averaging time is increased, the 
oscillations in Fig. 2 disappear. We note that the loca- 
tions of the boundaries between the different phases are 
not affected by the value of the velocity averaging time. 
At incommensurate fields, there is enough dispersion in 
the velocity of the moving interstitial vortices that the 
coherent velocity oscillations are no longer distinguish- 
able. 

As Fd increases, the net vortex velocity in the SB 
phase increases linearly until Fo = 0.365, where there 
is an abrupt increase in (14). Fig. 2(a,b) shows that this 
increase coincides with a jump in (Vy) to a zero aver- 
age, indicating that the vortices are moving only in the 
x-direction on average. In Fig. 3(c) we illustrate the dis- 
ordered vortex trajectories that occur in this phase at 
Fn = 0.42. The vortices are continually depinning and 
being repinned, and the order in the vortex lattice is 



lost, as shown in Fig. 3(d). We term this the random 
(R) phase. It resembles random dynamical phases that 
have previously been observed for vortices in square pin- 
ning arrays when the interstitial vortices begin to depin 
vortices from the pinning sites ,29J . Figure 2 shows that 
there are pronounced random fluctuations in (Vx) and 
(Vy) in phase R, and also that (Vx) does not increase 
linearly with Fo but has a curvature consistent with 
Vx = {Fd - F^f/"^, where Fj^ = 0.365 is the thresh- 
old value for the SB-R transition. In the SB phase, the 
number of moving vortices is constant and is equal to the 
number of interstitial vortices, while in the R phase the 
number of moving vortices increases with Fd . 

At Fd — 0.53, the system organizes into a one- 
dimensional flowing state where the vortex motion is 
locked along the pinning rows, as shown in Fig. 3(e,f) 
for Fd = 0.65. The onset of this phase also coincides 
with the decrease of fluctuations in (14) and the loss of 
fluctuations in I^, as shown in Fig. 4(c). For Fd > 0.53, 
all of the vortices are mobile and Fig. 2(a) illustrates 
that the (Vx) versus Fd curve becomes linear again. We 
term this the moving locked (ML) phase, since the vor- 
tex motion is effectively locked along the pinning sites. 
When the vortices are rapidly moving, the pinning sites 
have the same effect as a flashirig one-dimensional trough 
that channels the vortices [sl, [111 . The vortices assume 
a smectic structure in the ML phase, since different rows 
have different numbers of vortices, resulting in the forma- 
tion of aligned dislocations. The ML phase is essentially 
the same state found in square pinning arrays at high 
drives when B/B^ > 1.0 [33]. 

In previous studies of square pinning arrays with 
strong pinning, the initial motion of the vortices for 
B/Bij, > 1.0 occurred in the form of one-dimensional 
channels between the vortices trapped at the pinning 
sites 33|. In the honeycomb pinning array, similar flow 
occurs in the SB phase as shown in Fig. 3(a). For 
B/B^ < 1.5 in the honeycomb array, the initial inter- 
stitial flow for depinning in the x direction occurs via the 
flow of individual vortices in a zig-zag pattern around the 
pinned vortices. Since there is no dimer ordering for these 
fillings, no transverse response occurs for B/B^ < 1.5. 
For B/Bcf, > 1.5, the interstitial vortices begin to form 
dimer states when two interstitial vortices are captured 
in a single large interstitial site. The dimers can lower 
their orientational energy by aligning with each other in 
both the ground state and the moving states. Dimers 
can only remain aligned in the moving state if they are 
channeling along one of the symmetry axes of the pin- 
ning lattice. If the dimers were to move strictly in the 
x-direction, they would be forced directly into the pinned 
vortex in the pinning site to the right of each large in- 
terstitial site. This would destabilize the rodlike dimers. 
Instead, the dimers maintain their integrity by moving 
along ±30° to the x-axis. Within the moving state, if one 
of the dimers were to move along -|-30° while the remain- 
ing dimers were moving along —30°, the two interstitial 
vortices comprising the first dimer would be forced close 
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together, destabilizing the dimer state due to the repul- 
sive vortex-vortex interactions. Instead, all of the dimers 
move in the same direction. 

The SB-R transition occurs when the combined forces 
on the pinned vortices from the external drive and the 
moving dimers are strong enough to depin the pinned 
vortices. At the closest approach in the x direction be- 
tween a dimer and a pinned vortex, the frontmost dimer 
vortex is a distance ao/2 from the pinned vortex and 
the rear dimer vortex is a distance Soq from the pinned 
vortex, where Oo is the lattice constant of the undiluted 
triangular pinning lattice. In addition to the force from 
the dimerized vortices, the pinned vortex experiences an 
opposing force from the neighboring pinned vortex a dis- 
tance ao away. In a simple approximation, the driving 
force needed to depin a vortex at a pinning site is thus 
Fd ^ Fp- [{Kiiao/2) + ifi(3ao/2)) - ifi(ao)]. Set- 
ting Fp = 0.85 gives Fn = 0.41, close to the value of 
Fd = 0.37 for the SB-R transition in Fig. 2. Once the 
pinned vortices depin, the system enters the random (R) 
phase, and since Fd is still considerably less than Fp, it 
is possible for vortices to be pinned temporarily in phase 
R. 

Studies of square pinning arrays have shown that after 
the onset of a random dynamical phase, the vortices can 
organize into a more ordered phase of solitonlike pulse 
motion along the pinning rows, followed by a phase in 
which all of the vortices channel along the pinning rows 
[S^ . At the transition to the one-dimensional pulse like 
motion, a larger fraction of the vortices are pinned com- 
pared to the random phase, so a drop in (V^) with in- 
creasing Fu occurs, giving a negative differential con- 
ductivity. In the honeycomb pinning arrays for the pa- 
rameters we have chosen here, we do not observe one- 
dimensional pulse motion or negative differential con- 
ductivity for driving along the x-direction. For the one- 
dimensional pulse motion or the ML phase motion seen 
in Fig. 3(e,f) to occur, the vortices must be moving at 
a sufficiently high velocity for the pinning sites to act 
like a flashing trough. When the vortices move along 
the pinning rows, the vortex lattice structure adopts a 
highly anisotropic configuration which would be unsta- 
ble at Fu — 0. During the period of time when a vor- 
tex passes through a pinning site, the vortex is pulled 
toward the center of the pinning row, which stabilizes 
the one-dimensional motion. When the vortex is mov- 
ing between the pinning sites, it can drift away from the 
one-dimensional path until it encounters another pinning 
site. In Ref. [3^1, it was shown that for square pinning 
arrays, increasing the pinning radius Rp stabilized the 
one-dimensional flow down to lower values of Fd . In the 
honeycomb pinning array, the one-dimensional flow is less 
stable due to the fact that the vortices must move over 
the much wider large interstitial site, giving the vortices 
more time to drift away from the pinning row. Since this 
means that a larger value of Fjj is required to stabilize 
the one-dimensional motion, it should be more difficult 
in general to observe the onset of one-dimensional soli- 
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FIG. 5: The power spectra S{i') of the x component of the 
velocity Vx {t) for the three phases in Fig. 4. (a) The SB phase 
at Fd =0.25 shows a pronounced narrow band noise signa- 
ture, (b) The R phase at Fd = 0.42 has a broad band noise 
signature, (c) In the ML phase at Fd = 0.65, a number of 
different frequencies are present due to the fact that different 
rows of the vortices move at different velocities. 



tonlike motion or negative differential conductivity in the 
honeycomb pinning arrays than in the square pinning ar- 
rays. 



A. Fluctuations and Noise Characteristics 

In order to characterize the moving phases more quan- 
titatively, in Fig. 4 we show time traces 14 (t) and Vy (t) of 
the vortex velocities at fixed Fd for the different phases 
for the system in Fig. 2. In the symmetry broken (SB) 
phase at Fd = 0.25, shown in Fig. 4(a), Vx is greater than 
\Vy\ by tan(30°) or about 1.7. Here both components of 
the velocity show a pronounced oscillation which arises 
when the interstitial vortices move in a coherent fashion 
over the periodic potential substrate created by the im- 
mobile vortices in the pinning sites. In Fig. 5(a), we plot 
the corresponding power spectrum S{i') of Vx obtained 
from 
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There is a pronounced peak in S{v) at the frequency of 
the velocity oscillation in the SB phase, indicating that 
mode locking effects could appear at B/B,jj = 2.0 when 
the symmetry breaking flow occurs. In square pinning 
arrays, experiments (17i] and simulations j3^] revealed 
Shapiro step-like mode locking of interstitial moving vor- 
tices at B/Bfj, — 2.0. In the honeycomb lattice, since 
there is also a strong oscillation in Vy in the SB phase, we 
expect that transverse mode locking could occur if an ad- 
ditional ac drive is applied in the y-direction. Such mode 
locking would appear as steps in both and (Vy) ver- 
sus Fo in the SB phase. Transverse phase locking, which 
produces steps that are distinct from Shapiro steps, has 
been observed for the motion of vortices in square arrays 
[st! ]. In general, if the vortices already have an intrin- 
sic velocity oscillation in the transverse direction, then 
pronoimced transverse phase locking is possible. 

In Fig. 4(b) we plot the time trace of Vx and Vy for 
the random (R) phase at Fn — 0.42. In this case (\^) = 
0, and although both Vx and Vy show fluctuations, no 
oscillations or washboard frequencies appear. In Fig. 5(b) 
we show the corresponding S{v) for Vx^ where we find 
a broad band noise feature consistent with disordered 
plastic flow [l,[l,[ll|. Since there are no coherent velocity 
oscillations, mode locking should be absent in the random 
(R) phase. 

In Fig. 4(c) we plot Vx and Vy in the ML phase at 
Fo = 0.65, where Vx has been shifted down by a factor 
of 3 for clarity. There are no visible fluctuations in Vy due 
to the one-dimensional nature of the flow, but there are 
small periodic oscillations in Vx generated by the motion 
of the vortices over the periodic substrate. Due to the 
fact that different one-dimensional rows contain different 
numbers of vortices, producing dispersion in the vortex 
velocities, the oscillation in Vx is not as pronounced as 
in the SB phase. The corresponding power spectrum in 
Fig. 5(c) contains a rich variety of peaks due to the wide 
range of frequencies present in this phase. The main peak 
is smaller in magnitude than that found for the SB phase. 
As Fd increases, the frequency at which the first peak 
occurs also increases. It should be possible to generate 
phase locking in the ML phase; however, it would likely 
not be as pronounced as in the SB phase. These results 
suggest that noise fluctuations can be a useful technique 
for exploring the presence of different dynamical phases 
in periodic pinning arrays. 



B. Dynamical Symmetry Breaking in the Pinned 
Phase 

As previously noted, the ground state at B/B^ = 2.0 
is three-fold degenerate. When the dimers are aligned at 
either -1-30° or —30° to the x-axis in the ground state, the 
subsequent SB flow is aligned in the same direction as the 
ground state. It is also possible for the dimers to align in 
the y-direction, as shown in Fig. 6(a). At Fd = 0.09, the 
dimers and the vortices in the pinning sites are shifted 
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FIG. 6: Vortex positions (filled circles), pinning site locations 
(open circles), and trajectories (black lines) for a system with 
B/B^ = 2.0 and Fp = 0.85 which started in the y-aligned 
ground state, (a) The pinned phase at Fd = 0.09. Here the 
dimers and the vortices in the pinning sites have all shifted 
slightly to the right compared to the ground state due to the 
applied drive, (b) A rotational instability occurs at Fd ~ 
0.11, when the vortices move in a manner that allows the 
dimers to align along —30° to the a;-axis. There is also a 
small shift of the vortices in the pinning sites, (c) The pinned 
state at Fd ~ 0.12 where the dimers are aligned in the new 
—30° direction. 
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FIG. 7: The velocity (Vx) (upper curve) and {Vy) (lower 
curve) vs Fd for the system in Fig. 6. The rotational insta- 
bility seen in Fig. 6(b) appears as a positive peak in (Vx) and 
a negative peak in (Vy) just above Fd ~ 0.11. The system 
remains pinned until around Fd = 0.14. 



slightly to the right due to the applied drive. As Fu is 
further increased, a symmetry breaking transition occurs 
within the pinned phase. For F^ < 0.11 the dimers re- 
main aligned the y-direction; however, at Fo ~ 0.11, the 
rotational instability illustrated in Fig. 6(b) occurs. The 
dimers rotate in such a way that they end up aligned in 
the —30° direction. The interstitial vortex at the bottom 
of the dimer moves in the +x direction and by a smaller 
amount in the +y direction, while the vortex at the top 
of the dimer moves in the —y direction and by a smaller 
amount in the —x direction. There is also a slight shift 
of the vortices in the pinning sites that are closest to the 
bottom of each dimer. In Fig. 6(c) the rotation process 
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FIG. 8: Vortex positions (filled circles) and pinning site loca- 
tions (open circles) for the system in Fig. 2 at (a) B/B^ = 1.91 
and (b) B/B^ = 2.08. At these fillings, the long range orien- 
tational ordering of the dimer state is lost. 



is completed and the dimers are aligned in a new direc- 
tion, —30°. The vortices remain pinned until F]j = 0.14, 
at which point the system enters the SB phase. At finite 
temperatures, the dimer realignment occurs at even lower 
values of Fjj . The rearrangement can also be observed as 
a jump in (Vx) and (Vy) as shown in Fig. 7, where there 
is a positive spike in (Vx) and a negative spike in (Vy) 
near Fu — 0.11, in agreement with the motion shown in 
Fig. 6(b). We term this a dimer polarization effect since 
the driving force induces an alignment of the dimers. In 
runs with slightly different initial conditions, the dimers 
may align along the -1-30° direction with 50% probability. 



C. Dynamics for 1.5 < B/B^ < 2.5 

We next consider the effect of changing the vortex den- 
sity for fillings where interstitial dimers are present and 
the SB phase occurs. In Fig. 8(a) we illustrate the vortex 
positions for B/ = 1.91, where a mixture of monomers 
and dimers appear in the large interstitial sites. At this 
filling, the overall orientational ordering of the dimers is 
lost in the ground state, and the dimers are oriented only 
in local patches. For B/Bfp > 2.0, a mixture of intersti- 
tial dimers and trimers is present, as shown in Fig. 8(b) 
for B/Bcf, = 2.08, and the orientational ordering is again 
lost. In Ref. [5^ it was shown that the SB state still 
occurs at incommensurate fields as long as some dimers 
are present. If Fjj is suddenly increased from zero to a 
finite value at which only the interstitial vortices depin, 
the moving state for the incommensurate fields organizes 
into a dynamically symmetry broken state where all of 
the dimers flow along -f 30° or —30°. At the incommensu- 
rate fields, only the dimers undergo dynamical symmetry 
breaking; the monomers and trimers continue to move in 
the direction of the drive, with some fluctuations in the 
transverse direction. 

In Fig. 9 we plot {Vy) for the system in Fig. 2 at 
B/B^ = 1.89, 1.94, and 2.5. In Figs. 9(a,b), the same 
four phases described above are labeled. The SB phase 




FIG. 9: (Vy) versus Fd for (a) B/B^ = 1.89, (b) B/B^ = 
1.94, and (c) B/B^ = 2.5 for a system with the same param- 
eters as in Fig. 2. P: pinned phase; SB: symmetry broken 
phase; R: random phase; ML: moving locked phase. 



has opposite sign in Fig. 9(a) and Fig. 9(b); the dy- 
namical symmetry breaking can occur in either direction 
since there is no symmetry breaking in the ground state. 
If slightly different initial conditions are used, such as by 
changing the initial annealing procedure, the dynamical 
symmetry breaking has equal probability to occur in the 
positive or negative direction, as shown previously [58| . 
In Figs. 9(a) and (b) the initial portion of the SB phase 
has fluctuations in (Vy) due to the fact that we are in- 
creasing Fo at a finite rate and there is a transient time 
for the moving state to fully organize into the SB state, 
as studied previously ^1]. The transient time increases 
as \B/ B^ — 2.0\ increases. If we decrease AFd, the fluctu- 
ations at depinning are reduced; however, the boundary 
between the phases does not shift. 

In Fig. 9(c) at B/B^ = 2.5, (Vy) = since there are 
only trimer states present. The large oscillations in {Vy) 
occur when t he sy stem forms a completely ordered trimer 
ground state [50| and the vortex motion is highly coher- 
ent, similar to the effect seen in Fig. 2. For the rate at 
which we sample and average (Vy) versus F^, the peri- 
odic fluctuating vortex velocity is visible. For Fd > 0.5 
the system enters a partially moving locked phase where 
a portion of the vortices move along the pinning rows. 
There are, however, too many vortices to form straight 
one-dimensional chains of the type shown in Fig. 3(e) for 
B/Bip = 2.0. A buckling instability of the chains oc- 
curs since the amount of anisotropy that would occur if 
one-dimensional chains formed is too large for the vor- 
tex lattice to sustain. Instead, a partially moving locked 
(PML) phase forms with a disordered moving vortex lat- 
tice. This result is interesting since it indicates that mov- 
ing vortex phases do not always organize into ordered 
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FIG. 10: The dynamic phase diagram for Fd vs B/B^ high- 
lighting the diflterent dynamical phases. P: pinned, SB: sym- 
metry broken, R: random, ML: moving locked. Here Fp — 
0.85, Rp = 0.35A, and Up = 0.3125/A^ For B/B^ > 2.125, at 
high Fd the system forms a partially moving locked (PML) 
phase where not all of the vortices move along the pinning 
rows. 



FIG. 11: The critical depinning force Fc vs Fp for a system 
with B/B4, = 2.0, Rp = 0.35A, and Up = 0.3125/A^. For Fp < 
0.45 a partially pinned (PP) state forms which is illustrated 
in Fig. 12(a). For 0.45 < Fp < 1.75 the system forms the 
pinned (P) orientationally ordered dimer state such as that 
shown in Fig. 1(c). For Fp > 1.75, the pinned herringbone 
(PHB) state seen in Fig. 12(b) forms. 



states. A time trace of a PML state at fixed Fd shows 
much weaker velocity oscillations than those shown in 
Fig. 4(c) for the ML state. This suggests that phase 
locking with PML states will be very weak or absent. In 
previous work on phase locking for square arrays, it was 
shown that the phase locking is most pronounced at com- 
mensurate fields where the moving vortex structures are 
more ordered [13, HI] . 

By performing a series of simulations for varied B/B^, 
measuring the features in the velocity force curves and 
observing the vortex structures, we construct the dynam- 
ical phase diagram of Fjj vs B/B^ shown in Fig. 10. The 
depinning force marking the end of the pinned (P) phase 
shows peaks at B/B^ = 1.5, 2.0, and 2.5, corresponding 
to the commensurate and ordered ground states reported 
previously [HO]- The SB-R transition line is fairly flat 
as a function of B/B^ with an enhancement to higher 
values of Fo occurring near B/B^ = 2.0, while at the 
incommensurate fields, monomers or trimers create fluc- 
tuations that cause the vortices at the pinning sites to 
depin at slightly lower values of Fd. For B/B^ < 2.1, 
upon increasing Fo the random state organizes into a 
ML state where all the vortices move along the pinning 
rows as shown in Fig. 3(e), while for B/B^ > 2.1, the 
random state organizes into the PML state. The width 
of the random phase, as determined by the fluctuations 
in the velocity, increases and persists to higher values of 
Fd for increasing B/B^ at B/B^ > 2.1. For B/B^ < 1.5 
and B^ > 2.5, where dimers are no longer present, the 
SB phase is lost and a new set of dynamical phases arises 
which we discuss in more detail below. 
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FIG. 12: Vortex positions (filled circles) and pinning site 
locations (open circles) for the system in Fig. 11. (a) The 
partially pinned (PP) state at Fp = 0.25. The vortex lat- 
tice structure consists of a triangular lattice, and only half of 
the pinning sites are occupied, (b) The pinned herringbone 
(PHB) state at Fp — 2.0. Here the dimers do not all align in 
the same direction, but instead alternate in their alignment 
from row to row. 



D. Effect of Changing the Pinning Strength 

We next consider the effect of changing the pinning 
strength when B/B^ = 2.0. The four phases in Fig. 2 
occurred in a sample with i^p = 0.85. As i^p is varied, we 
find several different kinds of ordering within the pinned 
phase that affect the dynamics which occur at finite Fd . 
In Fig. 11 we plot the threshold depinning force Fc as 
a function of Fp. For Fp < 0.35 the pinning is weak 
enough that the vortex- vortex interactions dominate over 
the pinning energy and a nearly triangular vortex lattice 
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FIG. 13: Vortex positions (filled circles), pinning site loca- 
tions (open circles), and vortex trajectories for the system in 
Fig. 11. (a) Moving crystal (MC) state for Fp = 0.25 and 
Fd = 0.2. The partially pinned (PP) state from Fig. 12(a) 
depins elastically into a MC where half of the vortices move 
directly along the pinning rows and the other half of the vor- 
tices move in winding paths between the pinning sites, (b) 
Moving interstitial (MI) state for Fp = 2.0 and Fd = 0.225. 
The pinned herringbone (PHB) state depins into a MI state 
in which interstitial vortices move around the pinned vortices. 



forms, as shown in Fig. 12(a). In this arrangement, half 
of the pinning sites are still occupied, so the vortex lattice 
is partially pinned (PP) and there is a finite depinning 
threshold. This type of partially pinned vortex lattice 
was observed in previous simulations on honeycomb pin- 
ning lattices (soj , and similar partially pinned vortex lat- 
tice states have been predicted for square pinning arrays 
[s^ l and observed for metallic particles on periodic struc- 
tures [g^ . The depinning transition from the PP state is 
elastic, and all the vortices depin simultaneously to form 
the moving triangular crystal (MC) shown in Fig. 13(a). 
In the MC, half of the vortices move in one-dimensional 
paths along the pinning rows while the remaining vor- 
tices move through the interstitial regions with a small 
transverse oscillation. 

Figure 11 shows that the pinned ordered dimer state 
(P) forms for 0.45 < Fp < 1.75. Over the range 
0.45 < Fp < 0.55, the depinning from state P does not 
occur by the initial flow of the interstitial vortices into the 
SB phase, unlike the case shown earlier for Fp = 0.85. In- 
stead, for 0.45 < Fp < 0.55, both the interstitial vortices 
and the vortices at the pinning sites depin simultaneously 
and rearrange into the moving crystal (MC) state shown 
in Fig. 13(a). We also flnd a peak in Fc at Fp = 0.5. 
This peak occurs due to both the change in the pinning 
configuration and a change in the depinning process. For 
Fp < 0.45, only half of the pinning sites are occupied and 
the vortex lattice depins elastically. At 0.45 < Fp < 0.55, 
all of the pinning sites are now occupied in the P state, 
but the vortex lattice still depins elastically. The pinning 
energy that must be overcome to depin the lattice is in- 
creased compared to the PP state, leading to an increase 
in Fc- For 0.55 < Fp < 1.75, the depinning process is 
plastic and only the interstitial vortices flow at depinning 
to form the SB state. Since the plastic depinning process 



does not require pinned vortices to depin, the threshold 
force Fc drops, producing the peak in i^c at i^p = 0.5. 

For Fp > 0.7 in Fig. 11, the depinning threshold Fc 
slowly increases with increasing Fp, and a transition in 
the pinned vortex structure occurs at Fp — 1.75. For 
Fp > 1.75, the dimers in the pinned state are no longer 
aligned but form a pinned herringbone (PHB) type struc- 
ture such as that shown in Fig. 12(b), where the dimers 
are tilted in the same direction in one row and tilted 
in the opposite direction in the adjacent rows. Herring- 
bone ordering of dimers has previously been observed for 
colloidal dimers on triangular lattic es l 33l| and for vor- 
tices in kagome arrays at 5/3 filling [50|. At depinning, 
the PHB state does not form a SB phase but instead 
forms the winding interstitial phase shown in Fig. 13(b). 
The dimers break apart into two monomers, with one 
monomer passing around the pinned vortices in the pos- 
itive y direction and the other monomer passing the 
pinned vortices in the negative y direction. We term 
this state the moving interstitial (MI) phase. 

In contrast to the herringbone state, the aligned dimer 
or ferromagnetic ordering of the P state occurs when the 
pinned vortices adjacent to the dimers are able to un- 
dergo a periodic distortion within the pinning sites, re- 
ducing the interaction energy between the pinned and 
interstitial dimer vortices and permitting the dimer align- 
ment. If the dimers are aligned along -1-30°, as in Fig. 5(a) 
of Ref . , the two pinned vortices closest to each inter- 
stitial vortex in the dimer both splay outwards away from 
the -1-30° direction. When Fp is increased, the pinned 
vortices are pulled toward the center of each pinning site 
and are no longer able to distort in order to accommodate 
the aligned dimer state, so the pinned herringbone (PHB) 
state forms instead. This result suggests that there may 
be other types of ground state ordering for vortices in 
honeycomb and kagome arrays in addition to those that 
have been reported previously. It may be possible to use 
the size and shape of the pinning sites as a means of 
controlling the type of crystalline structure that forms. 

In Fig. 14, we plot (14) versus Fd for Fp = 0.25, 1.125, 
and 2.125. Figure 14(a) shows the elastic depinning pro- 
cess for the PP state at Fp = 0.25, which moves directly 
into the MC phase after depinning. At Fp = 1.125 in 
Fig. 14(b), (Vx) increases linearly with Fd through the 
SB phase. The slope of (14) increases in the R phase, 
and the velocity-force relationship becomes linear in the 
ML phase. For Fp = 2.125, Fig. 14(c) indicates that 
the depinning occurs in two steps. The first depinning 
transition of the interstitial vortices only takes the system 
from the P phase into the moving interstitial (MI) phase, 
while at the second depinning transition, the pinned vor- 
tices depin and the sample enters the ML phase. Unlike 
the behavior at Fp = 1.125 in Fig. 14(b), at Fp = 2.125 
the intermediate random (R) phase is lost and is replaced 
by a sharp jump into the ML phase. 

By conducting a series of simulations we construct the 
dynamical phase diagram as a function of Fp and Fd , as 
shown in Fig. 15. At high Fd, the MC phase forms for 
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FIG. 14: (14) vs Fd for B/S^ = 2.0, Rj, = 0.35A, and 
Up = 0.3125/A . (a) At Fp = 0.25, there is an elastic depin- 
ning transition between the partially pinned (PP) state and 
the moving crystal (MC) state, (b) At Fp = 1.125 the four 
dynamical phases arc present. P; pinned; SB: symmetry bro- 
ken; R: random; ML: moving locked, (c) At Fp = 2.125 the 
pinned herringbone (PHB) state depins into the moving in- 
terstitial (MI) state illustrated in Fig. 13(b). The transition 
between the MI and the moving locked (ML) state is much 
sharper than the R to ML transition shown in (b). 



Fp < 0.45, while for Fp > 0.45 the ML phase appears 
instead. The SB phase exists for 0.55 < Fp < 1.75, and 
the SB-R boundary shifts to higher Fd with increasing 
Fp until it terminates at Fp = 1.75. For Fp > 1.75, the 
PHB state occurs at low drive, and the system depins 
into the Ml phase. The MI phase also extends as far 
down as Fp = 1.5, where the system passes from the SB 
phase into a narrow window of the random (R) phase 
with increasing Fu before the vortices organize into the 
MI phase. As continues to increase, the vortices at 
the pinning sites depin and the system passes through a 
second narrow window of the R phase until the vortices 
organize into the ML phase. At high Fp, the R phase 
becomes vanishingly small and the system passes directly 
from the MI to the ML phase. The transition into the 
ML state increases linearly with increasing Fp, while the 
depinning force saturates with increasing Fp. 
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FIG. 15: The dynamic phase diagram of Fr> vs Fp for 
BjB^ = 2.0, Bp = 0.35A, and Up = 0.3125/A^. The three 
pinned phases arc: PP, the partially pinned phase shown in 
Fig. 12(a); P, the pinned ordered dimer phase in Fig. 1(c); and 
PHB, the pinned herringbone phase illustrated in Fig. 12(b). 
The dashed line separates the moving crystal phase (MC) 
shown in Fig. 13(a) from the moving locked (ML) phase of 
Fig. 3(c). The symmetry broken (SB) phase illustrated in 
Fig. 3(a) occurs at intermediate vales of Fp, while the mov- 
ing interstitial (MI) phase shown in Fig. 13(b) forms at higher 
values of Fp. The random (R) phase is illustrated in Fig. 3(c). 




FIG. 16: The dynamic phase diagram of Fd vs Rp for 
B/B^ = 2.0 and Fp = 0.85. PP: partially pinned phase; 
P: pinned dimer phase; MI2: moving interstitial phase 2; SB: 
symmetry broken phase; R: random phase; MC: moving crys- 
tal phase; ML: moving locked phase. For Rp > 0.2A we ob- 
serve the same phases illustrated in Fig. 15 at Fp = 0.85. The 
curves do not extend above Rp = 0.55A since for Rp > 0.55A, 
multiple vortex pinning at individual pinning sites occurs. For 
-Rp < 0.2A, a partially pinned phase appears and the initial 
depinning is into a new moving interstitial phase termed MI2, 
illustrated in Fig. 17(a). The upper dashed line separates the 
MC phase from the ML phase. 
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FIG. 17: (a) Vortex positions (filled circles), pinning site 
locations (open circles) and vortex trajectories (black lines) 
for the system in Fig. 16 showing the moving interstitial 2 
phase (MI2) at Rp = 0.15A and Fp = 0.2. (b) (K) vs Fd for 
the same system. A sharp transition from the MI2 phase to 
the random (R) phase occurs. 
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E. Changing Rp and 

We next examine the effects of changing the pinning 
radius in a system with fixed Fp — 0.85 and B/B^ — 2.0. 
In Fig. 16 we show the dynamic phase diagram for Fo 
versus Rp obtained from a series of simulations. For 
0.2A < Rp < 0.55A, we find the same three moving 
phases, SB, R, and ML, as in Fig. 2. For Rp > 0.55A, the 
pins are large enough to permit double vortex occupancy 
at the individual pinning sites; in this case, a new set 
of phases appears which we do not consider in this work. 
For Rp > 0.2A, the R-ML transition decreases in Fjy with 
increasing Rp since the larger pinning sites make it eas- 
ier for the vortices to localize along the pinning rows and 
flow in the one-dimensional motion of the ML phase. For 
Rp < 0.2A at low Fd, we find the partially pinned PP 
state illustrated in Fig. 12(a). The onset of the PP phase 
coincides with a drop in Fc at Rp « 0.2A. Unlike the PP 
phase that occurs at low Fp in Fig. 15, which depins elas- 
tically, the PP phase at small Rp depins into a moving 
interstitial phase that is distinct from the moving inter- 
stitial phase shown in Fig. 13(b). In Fig. 17(a) we illus- 
trate the vortex trajectories in the phase which we term 
the moving interstitial phase 2 (MI2). The vortices flow 
in winding interstitial channels; however, unlike the MI 
phase, in the MI2 phase only half of the pinning sites are 
occupied. In Fig. 17(b) we plot (Vx) vs Fu for a system 
with Rp = 0.15A. A clear two-step depinning transition 
occurs, with (Vx) increasing linearly with increasing Fd 
in the MI2 phase. At high Fd and Rp < 0.2A, the pin- 
ning sites are too small for the ML phase to occur, and 
instead the vortices flow in the MC phase illustrated in 
Fig. 13(a). 

We next consider samples with flxed B/B^ = 2.0, 
Rp = 0.35A, and Fp = 0.85, but vary the value of i?^ 
by changing the pinning density Hp. This alters the av- 
erage spacing between neighboring vortices. Up to this 
point we have used Up = 0.3125/A^. In Fig. 18 we illus- 
trate the dynamic phase diagram for Fd versus Up. As 
Up increases, Fc increases since the depinning of the in- 



FIG. 18: The dynamic phase diagram Fd versus pinning den- 
sity Hv, which determines B^, at B/B^ — 2.0, Fp — 0.85, and 
Rp — 0.35A. P: pinned dimer phase; MI: moving interstitial 
phase; SB: symmetry broken phase; R: random phase; ML: 
moving locked phase. The transition to the ML phase shifts 
to higher Fd with decreasing Up. The SB phase appears only 
for intermediate values of Up. 

terstitial vortices is determined by the potential created 
by the vortices located at the pinning sites, and as the 
vortex density increases, the depth of the interstitial pin- 
ning potential also increases. The R-ML transition shifts 
to higher Fd with decreasing Up. Since the distance be- 
tween the pinning sites increases with decreasing rip, the 
moving vortices spend less time in the pinning sites. This 
destabilizes the ML phase and the vortices must move at 
higher velocities for the effective trough potential to be 
able to stabilize the ML phase. For Up > 0.78/A^, the SB 
phase is lost, since at this pinning density the interactions 
between the interstitial vortices and the pinned vortices 
become sufficiently strong that the depinning of the in- 
terstitial vortices also causes the pinned vortices to depin. 
As a result, the system passes directly from the P phase 
to the R phase with increasing Fd- For Up < G.14/A^, the 
vortex-vortex interaction becomes weak enough that the 
system depins into the MI phase illustrated in Fig. 13(b). 
This also coincides with an increase in the value of Fd 
at which the R phase appears, since the moving intersti- 
tial vortices in the MI phase do not approach the pinned 
vortices as closely as they do in the SB phase. 

F. Effect of Finite Temperature 

We next consider the effect of finite temperature on 
the system in Fig. 2 with B/B^ = 2^Fp = 0.85, Rp = 
0.35A, and Up = 0.3125/A2. In Ref. we showed that 
a transition can occur at finite temperature in which the 
vortex n-mer states lose their orientational ordering and 
begin to rotate while remaining confined within the large 
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FIG. 19: The dynamic phase diagram of Fd vs thermal 
force F'^ for a system with B/B^ = 2.0, F^ = 0.85, np = 
0.3125/A^, and Rp = 0.35A. P: pinned dimer phase; SB: 
symmetry broken phase; R: random phase; and ML: moving 
locked phase. Inset: (14,) vs Fd at F^ = 1.25 where a smooth 
depinning transition occurs. 
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FIG. 20: Vortex positions (filled circles), pinning site loca- 
tions (open circles), and vortex trajectories (black lines) for 
the system in Fig. 19 at F'^ = 1.25. (a) At Fd = 0.1, there 
is appreciable creep of the interstitial vortices and the dimers 
have lost their orientational ordering and are rotating within 
the large interstitial sites, (b) At Fd = 0.2, the vortex flow 
is disordered and vortices are continually depinning and re- 
pinning, (c) At Fd = 1.0, all of the vortices are moving but 
there is still difi'usion in the direction transverse to the drive, 
(c) At Fd = 1.6 the vortices are about to enter the ML state 
where the motion is confined to one-dimensional channels. 



interstitial sites. This state was termed a vortex plastic 
crystal. In Ref. [Ill we demonstrated that the SB phase 
disappears in the vortex plastic crystal state. In Fig. 19 
we plot the dynamical phase diagram of Fd vs for 
the same system in Fig. 2. In our units, the dimers melt 
at = 1.0. Above the melting temperature, there is 
appreciable creep of the interstitial vortices as they hop 
from one large interstitial site to another, as illustrated 
in Fig. 20(a) for Fd = 0.1 and = 1.25. As Fd is 
further increased the system enters the random (R) phase 
shown in Fig. 20(b) for Fd = 0.25. At higher drives, the 
vortices begin to localize along the pinning rows in one- 
dimensional channels; however, there is still appreciable 
hopping from one row to another as shown in Fig. 20(c) 
for Fd = 1.0. At even higher drives, the ML is recovered 
as illustrated in Fig. 20(d) for Fd = 1.6. For F'^ > 
1.35, the ML phase is lost and the high Fd flow is in 
the random (R) phase shown in Fig. 20(c). In the inset 
to Fig. 20 we demonstrate that at — 1.25, the sharp 
features in the velocity force curve seen at F-^ = in 
Fig. 2 disappear. 



IV. DYNAMICS FOR DRIVING IN THE 
TRANSVERSE DIRECTION 

We now consider the case where Fd is applied along 
the ^/-direction, Fd = Fdy, for the same system as 
in Fig. 1(d) with Fp = 0.85, Rp = 0.35A, and Up = 
0.3125/A^. A different set of dynamic phases appear 
that are distinct from those found for driving in the x- 
direction. In particular, the SB phase is lost and the 
dimers align in the y-direction with the initial depinning 
occurring in one-dimensional interstitial flow paths. In 
Fig. 21 we plot the velocity force curves for B/ = 2.0, 
2.15, and 2.32. 

Figure 21(a) shows the three phases, pinned (P), 
one-dimensional moving interstitial (IDMI), and moving 
locked (ML) that occur at B/B^ — 2.0. In the P phase, 
if the ground state contains dimers which are aligned 
at either -1-30° or -30° to the polarization ef- 

fect is induced by the applied drive similar to the ef- 
fect discussed earlier. In this case, however, the dimers 
shift such that they are aligned in the y-direction. The 
one-dimensional moving interstitial (IDMI) state which 
appears above depinning is illustrated in Fig. 22(a) at 
Fd = 0.2, where the interstitial vortices move between 
the vortices in the pinning sites. Near Fd = 0.55 there 
is a sharp depinning transition for the vortices in the 
pinning sites. After this depinning transition occurs, the 
vortices very rapidly rearrange into a moving locked (ML) 
phase where the vortices move along the pinning sites, as 
shown in Fig. 22(b) for Fd = 0.8. In the ML phases 
for driving along the x and the y directions, the vortices 
travel in one-dimensional channels along a row or a col- 
umn of pinning sites, respectively. The pinning rows fol- 
lowed by the vortices for x-direction driving are evenly 
spaced in the y-direction, so the vortices flow through 
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FIG. 21: The velocity force curves (K,) vs Fo for Fd = Foy 
for the system in Fig. 2 with F-p = 0.85, i?p = 0.35A, 
and np — 0.3125/A^. P: pinned dimer phase; R: random 
phase; IDMI: one-dimensional moving interstitial phase, il- 
lustrated in Fig. 22(a); ML: moving locked phase, illustrated 
in Fig. 22(b). (a) BjB^ = 2.0. (b) B/B4, = 2.15, where 
an additional random flow phase occurs at depinning due to 
the presence of trimers. (c) B/B^ — 2.32, where there is no 
longer a IDMI phase. 
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FIG. 22: Vortex positions (filled circles), pinning site posi- 
tions (open circles), and vortex trajectories (black lines) for 
the system in Fig. 21(a) at B/B^ = 2.0. (a) One-dimensional 
moving interstitial (IDMI) phase at Fd ~ 0.2. (b) The mov- 
ing locked (ML) phase for Fo = Foy at Fd = 0.8. 




FIG. 23: The dynamical phase diagram for Fd vs B/B^ for 
Fd — FDy in a sample with Fp = 0.85, Rp = 0.35A, and 
Up — 0.3125/A^. P: pinned phase; IDMI: one-dimensional 
moving interstitial phase; R: random phase; ML: moving 
locked phase. For B/B^ > 2.3 the IDMI flow is lost. 



the centers of the pinning sites. In contrast, the pinning 
columns followed by the vortices for y-direction driving 
are unevenly spaced in the x direction due to the sym- 
metry of the honeycomb lattice. As a result, for the 
y-direction driving the vortices do not flow through the 
centers of the pinning sites, but are instead shifted to the 
right and left of the pinning sites in alternate columns, 
as seen in Fig. 22(b). This produces a more even spacing 
between the columns of moving vortices. Since different 
columns contain different numbers of vortices, the ML 
phase for y-direction driving has smectic type character- 
istics. The velocity-force curves for 1.5 < B/B^ < 2.0 
have the same general form as the curve in Fig. 21(a) and 
show the same three phases. 

For B/Bcf, > 2.0 the appearance of trimer states dis- 
rupts the IDMI flow since the trimers cannot align com- 
pletely in the y-direction. This produces random (R) 
vortex flow at depinning, as shown in Fig. 21(b), with 
diffusive vortex motion occurring along the a;-direction. 
There are more pronounced fluctuations in (Vy) in the 
R phase, and the velocity-force curve is nonlinear and 
lower than the extrapolated linear behavior in the IDMI 
phase that begins near Fu = 0.23. The trimers can block 
the one-dimensional channels of flow shown in Fig. 22(a), 
lowering the number of mobile vortices. At higher drives, 
the trimers depin, straighten into a linear configuration, 
and flow in the IDMI phase. At Fd « 0.5, the vortices 
in the pinning sites depin, resulting in a transition from 
the IDMI phase to the R phase. For sufficiently high 
drives, the ML phase forms. For driving along the x-axis 
at B/B^ > 2.0, we showed in Fig. 9(c) that the ML phase 
is lost due to a buckling transition of the one-dimensional 
chains of vortex motions, and that a partially moving 
locked (PML) phase forms instead when a portion of the 
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FIG. 25: The dynamic phase diagram Fd vs Fp for F_d = 
Foy at B/B4, = 2.0, Rp = 0.35A, and Up = 0.3125/ PP: 
partially pinned phase; P: pinned dimer phase; PHB: pinned 
herringbone phase; IDMI: one-dimensional moving intersti- 
tial phase; MC: moving crystal phase; ML: moving locked 
phase. A peak in the dcpinning threshold Fc occurs near 
Fp = 0.5 at the PP-P transition. At the P-PHB transition, 
Fc increases by a factor of three. The dashed line separates 
the MC phase from the ML phase. 



FIG. 24: (Vy) vs Fd for Fd = FdY for B/B^ = 2.0, 
Rp = 0.35A, and Up = 0.3125/A^ (a) At Fp = 0.35, there 
is a single step clastic dcpinning transition from the partially 
pinned (PP) phase to the moving crystal (MC) phase, (b) At 
Fp = 1.25 we find the pinned (P), one-dimensional moving 
interstitial (IDMI), and moving lattice (ML) phases, (c) At 
Fp = 2.25, there are sharp transitions between the pinned her- 
ringbone (PHB), one-dimensional moving interstitial (IDMI), 
and moving locked (ML) phases. 



vortices move through the interstitial regions. For driv- 
ing along the y-axis at B/B^ > 2.0, the ML state re- 
mains stable for much higher values of B/B^ than for 
the .T-axis case. A comparison of Fig. 22(b) and Fig. 3(c) 
shows that the interstitial region crossed by the vortices 
in the moving channels is not as wide for y-axis driving as 
for the a;-axis driving, resulting in more stable y-axis ML 
flow. As B/Bcf, is further increased, the random regime 
grows until the IDMI phase is completely lost, as shown 
in Fig. 21(c) for B/B^ = 2.35. In Fig. 23 we plot the dy- 
namical phase diagram for Fu versus B/B^, highlighting 
the onset of the different phases. The transition to the 
ML phase shifts to higher values of Fjj with increasing 
B/B^ since the ML vortex channels become increasingly 
anisotropic as the number of vortices in the sample in- 
creases. 



Dynamics as a function of Fp and Dimer 
Jamming 



We now consider the vortex dynamics in a system with 
fixed B/B4, = 2.0, Rp = 0.35A, and Up = 0.3125/A2 
for varying Fp with = ^dY- As noted above, for 
-Fp < 0.5 a partially pinned (PP) vortex lattice forms. In 
Fig. 24(a), the velocity- force curve for Fp = 0.35 shows 
that the depinning of the PP phase is elastic and occurs in 
a single step transition to a moving crystal (MC) phase, 
as seen earlier in Fig. 14(a) for driving along the x-axis. 
In the MC phase, half of the vortices move along the 
pinning sites. For 0.5 < Fp < 1.75, at low drive the sys- 
tem is in the pinned (P) phase of orientationally ordered 
dimers, and as the drive increases. Fig. 24(b) shows that 
the same IDMI and ML phases illustrated in Fig. 21(a) 
appear. The rapid rearrangement of the vortices from 
the IDMI phase to the ML phase results in a small jump 
near Fjj = 0.8 which marks the IDMl-ML transition. 
For strong pinning Fp > 1.75, Fig. 24(c) indicates that 
the ground state forms the pinned herringbone (PHB) 
phase illustrated in Fig. 12(b) . The same pinned state ap- 
pears for x-direction driving at strong pinning, as shown 
in Fig. 14(c). In Fig. 24(c), the velocity- force curve at 
Fp = 2.25 shows the abrupt nature of the depinning tran- 
sition from the PHB phase to the IDMI phase, which dif- 
fers from the smoother depinning transition that occurs 
from the P phase to the IDMI phase in Fig. 24(b). The 
depinning threshold increases markedly with increasing 
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FIG. 26: Vortex positions (filled circles) and pinning site 
locations (open circles) for the system in Fig. 24(c) at B / = 
2.0 and Fd ~ 0.56, just before depinning. Under the influence 
of the driving force which is applied in the y direction, the 
dimers align with the x direction and shift to the top of the 
large interstitial sites. Because the dimers are not aligned 
with the direction of the drive, a jamming phenomenon occurs 
which is responsible for the large increase in Fc seen in Fig. 25 
at Fp = 1.75. We call this the jammed state J. In (b) only 
the vortex positions are shown and it can more clearly be seen 
that the dimers are shifted in the positive y-direction. The 
vortex configuration in the jammed state is distinct from the 
pinned herringbone state. 



Fp once the system enters the PHB state. In Fig. 25 we 
plot the dynamic phase diagram for Fd versus Fp. Near 
the transition from the PP to the P phase, there is a peak 
in Fc similar to the peak observed at the PP-P transition 
for driving in the x-direction in Fig. 11. For Fp > 0.175 
the strong enhancement of the depinning threshold in the 
PHB state can be seen clearly. 

In Fig. 26(a) we illustrate the vortex positions just be- 
fore depinning for Fp = 2.25. Even though the drive is 
applied in the y-direction, the dimers have aligned with 
the a;-direction. When the dimers are oriented along 
the X-axis, they cannot fit through the easy-flow one- 
dimensional channel between the pinning sites, but in- 
stead are essentially jammed by the two pinned vortices 
at the top edge of the large interstitial site. In the or- 
dered dimer pinned (P) phase, the dimers all reorient in 
the same direction under an applied drive. In contrast, in 
the pinned herringbone (PHB) phase, the dimers rotate 
in opposite directions under an applied drive, so when 
the drive is applied along the y-direction the dimers end 
up aligning in the x-direction. In Fig. 26(b) only the 
vortex positions from Fig. 26(a) are shown to indicate 
more clearly the shift of the dimers in the positive y- 
direction. This vortex configuration, which we term the 
jammed (J) state, has a structure that is distinct from 
that of the pinned herringbone phase shown in Fig. 12(b). 
The jammed state configuration exists only in the pres- 
ence of the applied drive. For F^ ^ the dimers return 
to the pinned herringbone (PHB) state. In the jammed 
(J) state, the critical current is up to three times larger 
than in the state where the dimers are aligned in the y- 
direction. We also note that at incommensurate fields for 




FIG. 27: The critical depinning force in the a;-direction, F^f 
(open squares), and in the j/-direction, (filled circles), vs Fp 
for B/B4, = 2.0, Rp = 0.35A, and rip = Q.3125/A^ PP: par- 
tially pinned phase; P: pinned dimer phase; PHB: pinned her- 
ringbone phase; J: jammed state. In the PP phase, F^ = F^ , 
while in the P phase, F^ > F^ ■ A large enhancement of F^ 
occurs in the PHB phase when dimer jamming occurs. Inset: 
the ratio FJ^/J;^ vs Fp. The dashed line indicates F^/F^ = 1, 
where the depinning thresholds are equal. 



Fp > 1.75, the net vortex flow is reduced since some of 
the dimers align in the x-direction and effectively block 
the motion of other vortices along the j/-direction. 

In order to better characterize the enhancement of Fc 
in the jammed state, in Fig. 27 we plot the critical depin- 
ning force in the y-direction, F^ , and in the x-direction, 
F^, versus Fp. In the inset of Fig. 27 we show the ratio 
F^ /Fc versus Fp. In the partially pinned (PP) phase, 
F^ — Fc, while in the pinned aligned dimer (P) phase, 
F^ is slightly higher than Fy since the vortices can depin 
more readily into the IDMI phase in the y-direction. In 
the jammed state that forms from the PHB phase, F^ is 
3.1 times higher than for the same value of Fp. 



B. Effects of Changing Rp and B^ 

In Fig. 28 we plot the dynamical phase diagram Fd 
versus Rp for driving in the y-direction with B / = 2.0, 
Fp = 0.85, and Up = 0.3125/A2. For Rp > 0.2A, the sys- 
tem depins into the IDMI phase and makes a transition 
to the ML phase at higher drives. For Rp < 0.2A, the 
system forms the partially pinned (PP) phase where only 
half of the pinning sites are occupied. The PP phase de- 
pins into a moving interstitial phase (MI2Y) that resem- 
bles the MI2 state observed for driving in the x-direction 
in Fig. 17(b), where half the vortices depin while the 
other half remain pinned. The MI2Y phase is oriented 
90° from the MI2 phase. At Fd = 0.4 for Rp < 0.2A, 
the vortices at the pinning sites begin to depin and re- 
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FIG. 28: The dynamic phase diagram for Fd vs Rp with FIG. 29: The dynamical phase diagram for Fd vs np, which 

Fd = FdY, B/B4, = 2.0, Fp = 0.85, and Up = 0.3125/A^ 



PP: partially pinned phase; J: jammed state; P: pinned phase; 
MI2Y: y-direction moving interstitial phase 2; IDMI: one- 
dimensional moving interstitial phase; R: random phase; ML: 
moving locked phase. For Rp < 0.2A the system forms the PP 
phase. This depins into the MI2Y state for driving in the y- 
direction, which is similar to the MI2 state shown in Fig. 17(a) 
for driving in the x-direction. For 0.2A < Rp < 0.3A, the 
system forms the jammed J state shown in Fig. 26. 



pin, giving a regime of the random (R) phase until Fd 
becomes large enough for all the vortices to depin into 
the ML phase. For 0.2A < Rp < 0.3A, the jammed (J) 
state discussed in Fig. 26 occurs due to the formation 
of dimers aligned in the x-direction, which is associated 
vi^ith a marked increase in Fc. As discussed earlier, the 
pinned herringbone (PHB) phase and jammed (J) state 
occur when Fp becomes high enough that the vortices in 
the pinning sites cannot shift to allow for dimcr order- 
ing to occur. Similarly, as Rp is reduced, the vortices 
in the pinning sites have less room to adjust for dimer 
ordering, so the PHB state forms. The jamming also 
produces the counterintuitive effect that as Rp increases 
above Rp = 0.35A, the depinning threshold decreases. 



In Fig. 29 we show the dynamical phase diagram for Fjj 
versus Up, which determines the value of B^, for B/B^ = 
2.0, Rp = 0.35A, and Fp = 0.85. As Up increases, the 
critical depinning force into the IDMI phase increases 
since the repulsion from the pinned vortices experienced 
by the interstitial vortices increases as the average vortex- 
vortex spacing decreases. The transition from the IDMI 
phase to the ML phase shifts to higher values of Fd as 
Up decreases since the distance between the pinning sites 
which stabilize the ML flow increases. 



determines B^, for driving in the j/-direction with B/B^ 
2.0, Rp = 0.85A, and Fp = 0.85. P: pinned ordered dimer 
phase; IDMI: one-dimensional moving interstitial phase; ML: 
moving locked phase. 



V. DISCUSSION 

Our results are for honeycomb pinning arrays where 

it was shown in previous work that n-merization of the 
interstitial vortices into vortex molecular crystal states 
occurs for B/B^ > 1.5. Many of the dynamical effects 
presented in this work are due to the n-merization ef- 
fect. In kagome pinning arrays, similar types of vortex 
molecular crystal states appear, so we expect that many 
of the same types of dynamic phases described here will 
also occur for kagome pinning arrays, although we do ex- 
pect that there will be certain differences as well. In 
the kagome pinning array, the vortex dimer state ap- 
pears at B/B^ = 1.5 and has a herringbone ordering 
even for large, weak pins. There are no easy flow chan- 
nels along ±30° to the a:;-axis, so the symmetry breaking 
flows should be absent. Additionally, since there is no 
easy flow channel in the t/-direction, the anisotropic de- 
pinning dynamics may be different as well. 

We have only considered B/B^ < 2.5 in this work. 
At higher fields, a wide array of vortex molecular crystal 
states occur that should also have interesting dynamical 
phases. Since the low matching fields are more robust, 
observing the dynamics near these low fields experimen- 
tally is more feasible. Although our results are specifi- 
cally for pinning sites with single vortex occupation, sim- 
ilar dynamics should occur if the first few matching fields 
have multiple vortices at the pinning sites. In this case, 
the effective dimerization of the interstitial vortices would 
be shifted to higher magnetic fields. 

Although true phase transitions are associated only 
with equilibrium phenomena, the nonequilibrium phases 
considered here have many analogies to equilibrium 



18 



phases. For example, several of the transitions between 
the nonequilibrium phases have a continuous type be- 
havior, while in other cases the transitions are sharp, 
indicative of a first order nature. Future studies could 
explore the possible emergence of a growing correlation 
length near the transitions to see whether they exhibit 
the true power law behavior associated with continuous 
phase transitions or whether they show crossover behav- 
ior. For transitions that exhibit first order characteris- 
tics, it would be interesting to prepare a small patch of 
pinning sites with different characteristics that could act 
as a nucleation site for one of the phases in order to un- 
derstand whether there is a length scale analogous to a 
critical nucleus size. 

We also note that the dynamics we observe should 
be general to systems with similar geometries and re- 
pulsively interacting particles. For example, in colloidal 
systems, square pinning arrays with flat regions between 
the pinning sites (mufiin-tin potentials) have been fab- 
ricated, and in these systems the interstitial colloids are 
much more mobile than in washboard-type pinning po- 
tentials. Honeycomb pinning arrays could be created us- 
ing similar techniques for this type of system. 

VI. SUMMARY 

We have shown that vortices in honeycomb pinning ar- 
rays exhibit a rich variety of dynamical phases that are 
distinct from those found in triangular and square pin- 
ning arrays. The honeycomb pinning arrays allow for the 
appearance of n-mer type states that have orientational 
degrees of freedom. We specifically focused on the case 
where dimer states appear. At B/B^ = 2.0, the dimers 
can have a ferromagnetic type of ordering which is three- 
fold degenerate. At depinning, the dimers can flow in 
the direction in which they are aligned. For the case 
of driving along the x-axis, the dimers flow at ±30° to 
the applied drive, giving a transverse velocity response. 
At incommensurate fields where dimers are present, even 
though the orientational order is lost, the moving states 
can dynamically order into a broken symmetry state 
where the vortices flow with equal probability at either 
-1-30° or —30° to the a;-axis. As the driving in the x- 
direction increases, there is a depinning transition for the 



vortices in the pinning sites, and the transverse response 
is lost when the vortices either flow in a random phase or 
channel along the pinning sites. As a function of pinning 
force, we find other types of vortex lattice ordering at 
zero driving, including a partially pinned lattice and a 
herringbone ordering of the dimers. These other order- 
ings lead to new types of dynamical phases, including an 
elastic depinning for weak pinning where all the vortices 
depin simultaneously into a moving crystal phase, and 
an ordered interstitial flow in which the moving dimers 
break apart. The transitions between these flow phases 
appear as clear steps in the velocity force curves, and we 
have mapped the dynamical phase diagrams for various 
system parameters. We also showed that the different 
phases have distinct fluctuations and noise characteris- 
tics. When the temperature is high enough, the dimer 
states lose their orientational ordering and begin to ro- 
tate within the interstitial sites. This destroys the sym- 
metry breaking flow; however, the moving locked phase 
can still occur at high drives. 

The transition in the vortex ground state ordering as 
a function of pinning force causes the critical depinning 
force for driving in the x and y-directions to differ. When 
driving along the y-direction, the initial depinning oc- 
curs in the form of one-dimensional interstitial channels, 
and at high drives the vortices can form an anisotropic 
moving locked phase. We find a large enhancement of 
the depinning force in the y-direction associated with the 
pinned herringbone phase when the dimers align in the 
a::-direction and creates a jamming effect. The jammed 
state can enhance the critical depinning force by a factor 
of three, and can also arise for decreasing pinning size. 
We expect that many of the general features we observe 
will carry over to the higher matching fields in the honey- 
comb pinning arrays and in kagome arrays since ordered 
n-mers states occur for the kagome lattice as well. 
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